Quasiparticle Band Structure and Density Functional Theory: 
Single-Particle Excitations and Band Gaps in Lattice Models 
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We compare the quasiparticle band structure for a model insulator obtained from the fluctuation 
exchange approximation (FEA) with the eigenvalues of the corresponding density functional theory 
(DFT) and local density approximation (LDA). The discontinuity in the exchange-correlation po- 
tential for this model is small and the FEA and DFT band structures are in good agreement. In 
contrast to conventional wisdom, the LDA for this model overestimates the size of the band gap. 
We argue that this is a consequence of an FEA self-energy that is strongly frequency dependent, 
but essentially local. 
PACS numbers: 71.10.-w, 71.15.Mb, 71.10.Fd, 71.45. Gm 



I. INTRODUCTION 

Hohenberg and KohnEI showed that the total energy 
of the interacting electron gas in an external potential 
is given by the minimum of a functional of the electron 
density, n(r). This semipal result was followed a year 
later by the Kohn-ShamB 'trick' which reduces the orig- 
inal problem to an auxiliary problem involving fictitious 
non-interacting particles moving in a one-body effective 
potential composed of the external potential, the familiar 
Hartree potential, and an exchange-correlation potential, 
^a;cMi which is a universal nonlocal functional of the spa- 
tially inhomogeneous density. This innovation, together 
with local approximations for the exchange-correlation 
potential (LDA), underlies electronic structure calcula- 
tions for a wide variety of bulk materials, surfaces, and 
superlattices. Properties that can be obtained from to- 
tal energy calculatioris are generally found to be accurate 
within a few percent .H 

A by-product of the Kohn-Sham (KS) approach is a 
set of single-particle eigenvalues. Without rigorous ba- 
sis (even for the exact density functional theory), these 
are often taken to be single-particle excitation energies. 
This interpretation, together with the LDA, leads to rea- 
sonably good agreement with experimentally determined 
single-particle energies for many materials ,u but encoun- 
ters difficulties in numerous other materials ranging from 
commonplace semiconductors to exotic f-electron inter- 
metallics. Here we report on calculations that explore the 
relation between quasiparticle energies, eigenvalues ob- 
tained from density functional theory (DFT), and those 
from a corresponding LDA. 



Notable among those materials for which LDA eigen- 
values do not provide a correct account of single-particle 
excitation energies are seemingly ordinary insulators such 
as Si and diamond for which the LDA eigenvalues lead to 
sizeable underestimates of band gaps. A more dramatic 
discrepancy occurs for the strongly-correlated heavy- 
electron materials. For this class of rare-earth and ac- 
tinide compounds, effective masses obtained from elec- 
tronic structure calculations in the LDA differ by one 
to two orders of magnitude fcom those observed in de 
Haas van Alphen experimentsJj Among the ground states 
exhibited by the heavy-electron materials are insulators 
characterized by small gaps to low-energy charge exci- 
tations, which display experimental signatures of heavy- 
electron metals upon doping. In contrast to 'ordinary' 
insulators, LDA electronic structure calculations seem to 
overestimate the size |©f energy gaps in these 'Kondo' 
insulators as observedia for Ce3Bi4Pt3. The insulating 
parent compounds of the high-temperature supercon- 
ductors are also examples of materials for which the 
LDA does not adequately describe electronic correlatians; 
the LDA erroneously predicts that these are metals.Q In 
the absence of a rigorous theoretical connection between 
the eigenvalues of density functional theory and single- 
particle excitation energies, it has become a matter of de- 
bate whether improved approximation to the exchange- 
correlation potential can remove discrepancies between 
calculated eigenvalues and observed single-particle exci- 
tation energies. _ 

The only rigorous connectionD between KS eigenvalues 
and quasiparticle energies is that the eigenvalue of the 
highest occupied KS wave function is equal to the chemi- 



cal potential, /x. This resuJ: 
band gap of an insulator .tli 



I 



may be used to calculate the 



E,^E. 



N+l 
cbm 



vbm 7 



(1) 



where E^^^^^ is the KS eigenvalue of the appropriate con- 
duction band minimum calculated for the N + 1-particle 
system and -E^„ is the KS eigenvalue of the appropriate 
valence band maximum for the (insulating) iV-particle 
system. Here, any ambiguity in /x for a semiconductor is 
removed by taking appropriate limits for an M particle 
system as M — > A^ + 1 from above and M ^ N from 
low. The "band gap problem" with the KS eigenvalued 
appears if the gap is expressed soley in terms of the eigen- 
values of the A^-particle system E^ , 
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where AYxc is a spatially constant discontinuity of the 
exchange-correlation potential in going from an N parti- 
cle system with an exactly filled valence band to an iV + 1 
particle system with exactly one electron in the conduc- 
tion band. Because (for a macroscopic solid) the electron 
densities in these two cases differ by 0{1/N), AVxc is in 
a sense maximally nonlocal: it depends on the total elec- 
tron count (integrated over the entire system) and on the 
specific periodic structure of the solid. 

Insofar as AVxc is small for a particular insulator, the 
major source of error in a LDA calculation of its gap lies 
in the approximation for the remaining part of Vxc, which. 
is taken to be local, possibly with gradient corrections.^] 
Otherwise, AVxc rnust be explicitly included in an N 
particle LDA calculation for the gap. 

There may be no universal solution to these problems; 
AVxc may be large for one class of systems and small for 
others, and the physics underlying the contribution of 
AVxc may also differ from case to case. To address this 
issue, one needs quasiparticle energies, true DFT eigen- 
values, and LDA eigenvalues for the same system. Un- 
fortunately, one or another of these are known exactly in 
only a few special cases, and there is no even remotely re- 
alistic case for which all three are known exactly. Hence, 
one is always comparing approximate results, and a ma- 
jor concern is whether these comparisons are meaningful. 

An essentially exact LDA is known for 3D continuum 
systems, but in this case neither the true DFT eigen- 
values nor the exact quasiparticle energies are known 
with any pxeliability. For this case, Godby, Schliiter, 
and Shamllil studied semiconductors and tried to obtain 
meaningful comparisons by working consistently within 
the GW approximation, which was used to calculate the 
LDA exchange-correlation potential, quasiparticle ener- 
gies, and the true exchange-correlation potential for Si, 
GaAs, and AlAs. They found that the discontinuity in 
the exchange-correlation potential accounted for much of 
the underestimate of the band gap in the LDA. 

Other groups have worked with ID model Hamiltoni- 
ans with a small number of lattice sites, for which es- 
sentially exact quantum Monte Carlo calculations are 



possible .I13J13 Problems with this approach include the 
relevance of results from ID (withJjts well-known patholo- 
gies) to more realistic systems,li3 the lack of accurate 
quasiparticle energies for comparison, and the question 
of what constitutes an apipropriate LDA for compari- 
son with the exact results.E3 These references find that 
the discpntinuitypiii the exchange-correlation potential is 
small,E£l or largeJlj j— . 

In the spirit of Godby, Schliiter and Sham,llll we 
have studied the DFT, LDA, and quasiparticle ener- 
gies for a 2D two-band model insulator on a bipar- 
tite lattice with Hubbard on-site interactions. In place 
of the GW approximation we use the fluctuation ex- 
change approximationta (EEA), a fully self-consistent 
conserving approximation,EJ to calculate quasiparticle 
band structures and the thermodynamic potential. Be- 
cause we work at finite temperatjures, we use Mermin's 
finite-temperature generalizationllS of the Hohenberg- 
Kohn theorems to provide the connection between true 
DFT and the FEA: we define our "fuU DFT" to repro- 
duce the site densities of the FEA and we define our LDA 
for the exchange-correlation potential from the FEA for 
the single-band Hubbard model, which plays a role anal- 
ogous to that of the homogeneous interacting electron 
gas in continuum DFT. We consider model parameters 
for which no fluctuation channel is strongly dominant, 
but for which spin fluctuations and fluctuations in the 
particle-particle channel provide larger contributions to 
the thermodynamic potential than do density fluctua- 
tions. 

Calculations for our model insulator show that the 
FEA band gap is slightly smaller than that for the DFT 
and that the band gap for the LDA is signiflcantly larger 
than the FEA gap. These relationships among gaps dif- 
fer signiflcantly from the conventional electronic struc- 
ture wisdom. We believe that the relationships we find 
are a consequence of a strongly frequency-dependent but 
essentially local self-energy. Our self-energies are calcu- 
lated using a propagator-renormalized theory in which 
the bare parameters of the Hamiltonijaii have been elimi- 
nated from thermodynamic potentialJlj This propagator- 
functional formalism casts the thermodyanmic potential 
as a functional of observable spectral functions .E3 From 
this viewpoint, the form of the self-energy is a spectral 
property with generalizable physical consequences that 
transcend the specific underlying Hamiltonian. The ubiq- 
uitous quasiparticles of Landau's Fermi liquid theory pro- 
vide a familiar example. They arise not only in liquid 
■^He and in metals, but appear also in self-consistent cal- 
culations based on Hubbard Hamiltonians in three di- 
mensions. In this spirit, we belive that the effects of the 
locality and non-locality of the self-energy transcend the 
underlying Hamiltonians from which they were deduced, 
and so the essential physics distilled from our results may 
be applied to the 'Kondo insulators,' such as Ce3Bi4Pt3, 
for which electronic structure calculations in the LDA 
overestimate band gaps.Q 

In the next section, we introduce single-band and 



bipartite-lattice Hubbard models and the fluctuation ex- 
change approximation. In Section III the construction 
of a suitable density functional theory and local density 
approximation is described. In Section IV, the LDA and 
FEA are compared. Our conclusions are presented in 
Section V. 



II. THE MODEL INSULATOR 

To explore the relation of the eigenvalues of density 
functional theory to single-particle excitation energies, 
we take a Hubbard model on a bipartite lattice for our 
model insulator. The lattice consists of two interpene- 
trating square (Bravais) lattices, 'A' and 'B;' an A-site's 
nearest neighbors are B-sites and its next nearest neigh- 
bors are A-sites. The Hamiltonian for our model is 
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where Cjcr{Cj„) annihilates (creates) a particle on site j, 
tij are hopping matrix elements, Vi is an on-site single 
particle potential that takes the value va {vb) for a site 
on the A (B) sublattice, and U is the intra-site repulsion 
experienced by two electrons on the same lattice site. 
We work in 2D and include nearest-neighbor t and next- 
nearest-neighbor txy hopping amplitudes to avoid known 
pathologies that result from Fermi surface nesting and 
van Hove singularities at the Fermi surface. 

We write the one-body potential for the A (B) sites 
in the suggestive form w^ = vq — Aq [vb = wo + Aq). 
This potential doubles the unit cell leading to a two- 
dimensional analog of the Sru-,Schrieffer, and Heeger 
Hamiltonian for polyacetylene.cll For [/ — 0, the single- 
particle excitation energies are 
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where 



7k = -2 t^y (cos k^ + cos ky), 
a^ = -At cosfcj;COsfcye'*^''"+'''''^/^. 



(5) 



For nonzero U , the exact solution of this model is un- 
known. We use the fluctuation exchange approximation 
(FEA), based on the propagator renormalized perturba- 
tion theory of Luttinger and Ward, to provide quasipar- 
ticle energies and the basis for constructing full and LDA 
density functional theories^ 

The Luttinger and WardE3 formula expresses the grand 
thermodynamic potential f2 as a functional of the fully 
renormalized Green's function G and the self energy E, 



Here 'Tr' indicates a sum over momentum, frequency, and 
sublattice argumentsjHj and Go is the Green's function 
for the noninteracting system ([/ = 0) which contains 
the bare band structure. When viewed as a functional 
of G and E, fJ is stationary with respect to independent 
variations around the true G and S for a given Go and 
[/, leading to the 'skeleton diagram' expansion for the 
self-energy, 



^^^'""^"2 JG(k,£„)' 
and to Dyson's equation, 

G-i(k,e„) = Go'(k,£„)-E(k,e„) 



(7) 



(8) 



The self-consistent solutions of Eqs. (0) and (H) com- 
pletely specify G and E, and together with Eq. (0) 
provide a self-consistent description of single-particle ex- 
citation spectra and thermodynamic properties. The 
FEA entails a specific prescription for ^[G] that includes 
contributions from the exchange of spin, density, and 
Cooper-pair fluctuations in addition to self-consistent 
contributions at first (Hartree Fock) and second order 
in the interaction. The propagator for our model is a 
2x2 matrix labeled by sublattice indices; a more de- 
tailed discussion of the FEA for the Hubbard model on 
a bipartite lattice, together with explicit expressions for 
^[G], appears in the Appendix. 

We have solved Eqs. (0) and (g) self-consistently on 
a 64 X 64_mesh using a parallel algorithm described 
elsewhere.Ej This is not our most accurate method for 
solving the equations of the FEA. The more accurate 
algorithm of Ref. |2^, which does not contain a high- 
frequency cutoff in the traditional sense, is essential 
for calculating temperature derivatives of the thermody- 
namic potential and for performing calculations near an 
instability to an ordered state. Since none of these con- 
ditions apply here, the algorithm of Ref. |^ suffices for 
our purposes. 

The A- and B-site densities needed for the density func- 
tional theory calculations described below are obtained 
directly from the diagonal components of the propagator. 



2T 



/ ^ '-'ii('^: £riji 



(9) 



where T is the temperature and N is the number of k- 
points. Because the FEA is a conserving approximation, 
this is equivalent to calculating the site densities from 
the thermodynamic potential. 



dn[T,VA,VB] 

dvi 



(10) 



Explicit calculations show that this internal self- 
consistency is achieved to better than 1%. 



n{T,^l) = -2^T [SG + ln(- G^^ +E)] + $[G]. (6) 



III. DENSITY FUNCTIONAL THEORY FOR 
LATTICE MODELS 

Following Mermin,E3 it is straightforward to show that 
the thermodynamic potential of our model is a stationary 
function of the A- and B-site densities. The thermody- 
namic potential (with arguments T and ^ suppressed) 
is 

"[{^j}] - 0[{"j}] + ^H[{nj}] + n.,,[{n,}] + r!e,t[{n,}], 

(11) 

where ^ext is the contribution from the external poten- 
tial, O -|- ^ext is the grand thermodynamic potential of 
a system of fictitious non-interacting particles with site 
densities {uj}, and fi/f, and U,xc are the Hartree and 
exchange-correlation potential contributions. For our 
model, the Hartree contribution to fi is 



VLH[nA,nB\ = - (n\^n\) [/, 



(12) 



and the contribution from the external potential Vj is 

'^ext{{nj\\=^Vjny (13) 



The condition that the grand thermodynamic potential 
be stationary with respect to variations in the density 
leads to the set of coupled equations. 



de[{nj}] d£L^xt\{njY\ d£Lii\{nj\\ dVLxc\{njY\ 
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A. Full Density Functional Theory 

The KS formulation at finite temperature invokes an 
auxiliary system of noninteracting particles in a one-body 
effective potential T4ff , that is defined by writing Eq. (|lj) 
in the form, 
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The thermodynamic potential for the auxiliary problem 

is 
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where the single-particle eigenvalues Cq. satisfy the 
Schrodinger equation 

(i^ + Feff[{n,}])Va=e„Va, (17) 

and K is the hopping Hamiltonian (first term of Eq. (^). 
The site densities at finite temperature are related to the 
KS wave functions through 



2E 
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where /(e) is the usual Fermi function. Because Eq. ( [T5| ) 
is also the KS variational condition for the auxiliary prob- 
lem with Veff regarded as an external potential, densities 
that result from the self-consistent solution of Eqs. (|l7| ) 
and (nS) are guaranteed to be identical to those of the 
original interacting system. 

For our model, V^j^cij^-j}] is unknown; we use the FEA 
site densities together with the KS formulation to find the 
exchange correlation potential, Vxc = {''^xct'^xc)- To do 
this, we find one body potentials v^f and w^^, for which 
the site-densities of the finite-temperature auxiliary KS 
problem are the same as the FEA site-densities, and in- 
vert the definition of VcS to obtain Vxc- A convenient 
expression (equivalent to Eq. (18[)) for the site-densities 
of the non-interacting problem is obtained by taking the 
trace of the site-diagonal components of the propagator 
for the non-interacting system given in Eq. (A2). The 
full DFT band structure can be extracted from Eqs. (H) 
using Vi = vlg. 

Because the FEA is a conserving approximation, the 
FEA grand thermodynamic potential has the same sta- 
tionary properties with respect to variations of G and S 
as the exact grand thermodynamic potential. As a conse- 
quence, the exchange-correlation potential that we have 
calculated is identical to that calculated from Sham's in- 
tegral equationEa for 14c- 



B. The Local Density Approximation 

For our model, the single-band Hubbard model plays 
the role analogous to that of the uniform interacting elec- 
tron gas in the formulation of the continuum LDA for 
electronic structure theory. So, we use Eq. (|ll|) to cal- 
culate ria:cM in the FEA for the single-band Hubbard 
model (the Hamiltonian of Eq. g with Vi = 0) with the 
same "Coulomb potential," U. To avoid possible confu- 
sion, we emphasize that in this case the electron density is 
the same at every site, and fixe becomes a simple function 
of the uniform density, n. This entails evaluating Eq. (g) 
with the self-consistently calculated FEA propagator for 
the single-band Hubbard model, evaluating the Hartree 
contribution to the thermodynamic potential, and rear- 
ranging Eq. ^n\ ) with ^ext = to extract ^xc- The LDA 
exchange-correlation potential is 
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(19) 



In this way the LDA exchange-correlation potential cor- 
responding to the full DFT is calculated in an approxi- 
mation equivalent to that of the FEA calculation. 

Given Vx^J^'^in), Eqs. (|l^ and (^ can be solved 
directly with external potentials vq and Aq, subject to 
the constraint that the total density agrees with that of 
the FEA calculation. As in the case of the DFT, Eqs. 



(^ with the potentials vj^^g determine the LDA band 
structure. 



IV. RESULTS 



We focus on single-band and two-band models with 
next-nearest-neighbor hopping matrix element t^y = 0.35 
and [/ = 3. This set of parameters avoids possible com- 
plications of dominant fluctuation contributions in any 
one specific channel as might arise rfts a consequence of, 
for example, Fermi surface nesting.EH For the model in- 
sulator (see below), the contributions to the grand ther- 
modynamic potential from the Hartree and second-order 
diagrams are the largest. Somewhat smaller spin, den- 
sity, and Cooper pair fluctuation contributions are all 
roughly the same order. For T = 0.02, explicit val- 
ues for the $-functionals in their respective channels are: 
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and $;{ = 
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3.31 X 10" 
The exchange-correlation contribution to the thermo- 
dynamic potential ^Ixc as a function of density for a 
single-band Hubbard model with next-nearest-neighbor 
hopping matrix element txy = 0.35, ^ = 3 and T = 0.08 
is shown in Fig. ^ (we henceforth measure energies in 
units of the nearest-neighbor hopping matrix element t). 
The solid line is a polynomial fit of the form. 




FIG. 1. The exchange-correlation contribution to the 
thermodynamic potential Q.xc as a function of density for the 
single-band Hubbard model in the fluctuation exchange ap- 
proximation (FEA) for f/ = 3 at two temperatures. The solid 
line is a fit to 77 points (see text) for T = 0.08. The widely 
spaced points are calculated for T = 0.02. 
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through 77 data points. The coefficients, presented in 
Table 1, lead to a fit with an error of less than 1% over 
the entire density range; the worst agreement occurring 
at low density. It is interesting to note, in contrast to 
the exchange-correlation energy used in by Ref. O, that 



TABLE I. Values of the coefiicients a^ in Eq. (Q for 
the exchange-correlation contribution to the grand thermo- 
dynamic potential as a function of density. These parameters 



are for [/ = 3, T = 0.08, and t 



xy 



0.35. 



02 = -1.32532 
as = -13.7277 
as = 21.9011 



as = -0.300919 
ae = 29.0402 
a<j = -7.11752 



a4 = 3.27475 
ar = -34.1589 
aio = 0.917821 



there-is no evidence for a significant n^^^ contribution to 
rixc.tJ Fig. ^ also shows flxdn) for T = 0.02 at several 
densities, from which we conclude that flxc{n) is weakly 
temperature dependent over the range 0.02 < T < 0.08. 
The function ^xciji) is to a reasonable approximation 
linear in U (see Fig. ||); the largest deviations occur 
at high density as might be expected. We calculate the 
local density approximation of the exchange correlation 
potential from Eq. (EQ). 




FIG. 2. The exchange-correlation contribution to the ther- 
modynamic potential ilxcin) for U — 3 (solid line; see previ- 
ous figure and text), (7 = 1 (■), and the U — 1 result scaled 
by a factor of 3 (a). 

We now compare the FEA single-particle excitations 
of an insulator to the eigenvalues obtained from the full 
DFT and those from the LDA. To this end we calcu- 
late the, fully self-consistent self-energy in the FEA for 
a bareca staggered potential Aq = 2.532. The quasipar- 
ticle band structure is obtained from the zeros of the 
real part of the denominator of the retarded Green's 
function G^. The fully renormalized retarded propaga- 
tor was constructed using Dyson's equation and the re- 
tarded self-energy. The latter was obtained from an ana- 
lytic continuation of the self-energy from the imaginary- 
frequency axis to the_real-frequency axis using N-point 
Fade Approximants.cj The resulting band structure for 
Ao = 2.532 and n = 1.000 (with n^^^ = 1.688 and 
Ug^^ = 0.312), shown in Fig. ^ for k along the high- 
symmetry path T—X—M—T in the square zone, is that of 
an insulator with an indirect gap Eg = E'^—EJj = 1.318. 
The corresponding full DFT band structure is found by 
evaluating the (bare) dispersion relations in Eq. (0) using 
on-site one body potentials v, A for which the site den- 



sities UA and n;, are equal to n^^^ and n^^^. Starting 
from the (renormalized) on-site single-body potentials of 
the FEA (w^^^, A^^^), a Newton-Raphson method al- 
lowed us to determine (to 1 part in 10^) v^^'^ = —0.745 
and A-°^^ = 1.368. The DFT band structure appears 
in comparison with that for the FEA in Fig. g. The 
fact that these band structures are essentially indistin- 
guishable suggests that the discontinuity in the exchange- 
correlation potential is small. In Fig. ffl we show the 
'exact' V^^'^'^(n) obtained from v^^^ and A-^^-^ for a 
range of densities about n = \ and observe that, as 
expected, V^^ and V^^ each show a small discontinuity 
|AT4c| = 0.026. Note that the small structure appearing 
at the high-density side [n — 1.00+) of the discontinuity 
is not physically significant, but is rather a numerical ar- 
tifact reflecting a finite high-frequency cutoff in the FEA 
calculation; this feature is systematically reduced as the 
number of Matsubara frequencies kept in the calculation 
is increased. E3 




FIG. 3. The band structure for k along high symmetry di- 
rections in the Brillouin zone from the FEA calculation (solid) 
compared with that from the DFT (o) for T = 0.08, U = Z, 
and txy — 0.35. 

To calculate the band structure in the LDA, we take 
yLDA fj.Q^-^ gqg^ Q ^^^ Q^ ^^^ gQ^yg gqg^ Q ^^^ 

(Eq) using a Newton-Raphson method to find on-site one- 
body potentials (yLDA^j^LDA-^ ^-^aX yield the same total 
density with the same bare staggered potential Aq as in 
the FEA calculation. The LDA band structure obtained 
from Eqs. (13) with {yi^DA^^LDA^ jg shown compared 
with that for the FEA in Fig. ||. For n « 1, dn/d^i is 
vanishingly small and a wide range of v'"^^ lead to a 
Aq and n equal to those of the FEA calculation. On the 
other hand, A^^ is determined to high accuracy by our 
procedure and is largely insensitive to changes in v^^ 
for n '--^ 1. Since T is much smaller than the (indirect) 
band gap resulting from A^^^, we take v^^^ so that 
the chemical potential lies in the middle of the gap as 
expected on physical grounds. We have explicitly verified 
that this choice leads to n — n^^'^ and Aq = A^-^^ to 
within the accuracy to which these quantities are known. 
The LDA and FEA band structures differ significantly; 
the LDA band structure cannot be corrected by uni- 
formly and independently shifting the bands by a con- 



stant amount. The LDA band structure cannot be cor- 
rected by using the "scissors operator." An examination 
of the spatial dependence of the FEA self-energy shows 
that the E^ are essentially local. The self-energy is es- 
sentially uniform in k and strongly frequency dependent. 
The excitation energies Ek are the poles of the retarded 
propagator, G^ . By combining the appropriate analytic 
continuation of Eq. (A2) with Dyson's equation, we ex- 
press G^ in terms of the components of the self-energy 
and find that the poles of G^ satisfy 



[Ek + Ao + /^ + 7k - SsB(£^k)] - 

[Qk + ^AsiEk)] K + SBA(Ek)] = 0. 



(21) 



Here a^ and 7k are defined as in Eq. (pi) and we have used 
Sy (k, e) '^ Sij(e). The i?k that solves Eq. ( pi| ) has a k- 
dependence that not only reflects of that of 7k and ak but 
also includes implicit additional k-dependence through 
the frequency dependence of the S^ . 
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FIG. 4. The 'exact' exchange-correlation potential V^c 
(see text) for A-(») and B-sites(A). The discontinuity is 
~ 0.02 and occurs at half fiUing. 

As is plainly evident from Fig. Q the LDA over- 
estimates the the band gap by about a factor of 3. 
In the light of LDA calculations for insulators, the re- 
sult Eg^^ > Eg^^ is perhaps surprising. Since the 
FEA and the self-consistent GW approximation can be 
viewed as propagator- functional theories, physical prop- 
erties may be viessied as a consequence of the form of 
spectral functionsEj rather than of the form of a particu- 
lar Hamiltonian. Considered in this way, results of con- 
tinuum calculations are relevant insofar as they lead to 
trends in self-energies and hence spectral functions. From 
this perspective, we believe that the explanation for ops. 
results lies in the work of Godby, Schliiter, and Sham.til 
They observed that, for the GW approximation, the non- 
locality of the self-energy tends to a widen the gap while 
frequency dependence tends to narrow the gap. Given 
the importance of the frequency dependence of S^ in 
determining the FEA band structure and given that the 
DFT band structure is obtained from dispersion relations 



for a non-interacting system, albeit with largely 'renor- 
malized' on-site potentials, the observed close agreement 
between these band structures may be of a rather subtle 
physical origin. This agreement may also be in part a 
reflection of a moderate interaction strength {U ~ 3/8 
of the bare, Aq — 0, bandwidth) and may degrade with 
larger U. It is also evident from the DFT calculation (see 
Fig. ^) that V^ and V^^ are of opposite sign for the den- 
sity range in which we are working. On the other hand, 
fixe of Fig. d leads to a V^J^^ < over the entire density 
range. To compensate, the LDA must increase the size 
of A'^^'^, which increases the gap. 
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FIG. 5. A comparison of the FEA band structure (solid) 
with that of the LDA (o), where the LDA is adjusted to yield 
the same bare staggered field and total density of the FEA 
calculation. 



V. CONCLUSIONS 

We have calculated the renormalized band structure 
and site densities for a model insulator in the fluctuation 
exchange approximation (FEA) and in corresponding 
full and approximate (LDA) density functional theories. 
The band structures calculated from the full density- 
functional theory are in good agreement with the 'true' 
band structures obtained from the poles of the retarded 
FEA propagator. The gap obtained from the DFT is 
slightly smaller than that of the FEA; the difference is 
accounted for by a discontinuity in V^c that is small in 
magnitude. It is also of opposite sign to that observed 
in GW calculations for 'real' ordinary insulators. We 
find that gap in the LDA is larger than that of the full 
DFT. The relationship among the gaps differs from that 
observed in continuum GW calculations for ordinary in- 
sulators. An examination of the FEA self-energy shows 
it to be essentially local, suggesting that this is a con- 
sequence of the frequency dependence of the self-energy. 
We speculate that such a relationship ought to be ob- 
served in the 'Kondo' insulators where it is expected that, 
owing to strong electronic correlations, the frequency de- 
pendence of the self-energy is more important than its 



spatial nonlocality. 
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APPENDIX: THE FEA FOR A HUBBARD 
MODEL ON A BIPARTITE LATTICE 

The single-particle propagator is central to the calcula- 
tion of the density and single-particle excitation spectra. 
It is convenient to adopt a matrix notation for the prop- 
agator in imaginary time for the Hamiltonian in Eq. (0), 

e,,(r, -r;,r) = - < r{^(r„T)^t(r^.^ q)} > , (Al) 

where T is the usual Wick's time ordering operator and 
Tpiri) is an annihilation operator acting at r^, where i is a 
sublattice index {i ~ A,B). Because of the translational 
invariance of the full Bravais lattice, the propagator is 
a function of a coordinate difference. Unlike Gaa, Gab 
is defined on a Bravais lattice that does not include the 
origin. We find it convenient for numerical calculations 
on parallel computers to introduce a propagator Gij de- 
fined so that Gu = Q.U, Gab{t!^,t) = Qab{y - c,t), and 
GsAiTc, t) = Qba{t^ + c,t), where c is a vector that joins 
a point on the A sublattice to one on the B sublattice. 
By this definition, the arguments of all components of 
the propagator range over the same (sub)lattice and the 
usual Fourier transform connects (r, r) space to (k, £„) 
for all components of Gij . 

Dyson's equation (Eq. (ph) should now be inter- 
preted as an equation for the (2 x 2) matrices G, Gq 
and E. Keeping only nearest-neighbor and next-nearest- 
neighbor hopping amplitudes, the propagator for the 
nonintcracting system in (k, e„) space is, 



Go\K£n) = 



ISr, 



7k- 



IC-f-j 



-ttk 

7k- 



Ao 



(A2) 



where "fk — Jk + v — /i, and ji^ and a^ defined as in Eq. 
d) above. „ 

The FEAlla takes $[{Gij}] to be the sum of the 
Hartree-Fock diagram, the single second-order diagram, 
and particle-hole and particle-particle bubble chains that 
describe exchanged density, spin-density, and (singlet) 
pair fluctuations. For the bipartite lattice, explicit ex- 
pressions for 'I' [{Gij}] are 
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2 ^^ 



vh n2 



vh x2 



ix^AV + {x'bbY + 2x7bX'ba (A3) 
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Ki = iTr{ln((l + xfA)i^ + x'bb) ^ x'abX'ba) 
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'Xaa ~ Xbb 
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+X^AA + Xbb 
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p/^ ph 



%, = iTr{ln((l + xT4)(l + xg'B) - xJbX^/^) 



pp pp 

~Xaa ~ Xbb 
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(A4) 



(A5) 



(A6) 



Here x^- and x^f S'l'^ particlc-holc and particle-particle 
susceptibility bubbles connecting lattice sites i and j and 
are related to the fully rcnormalized Green's function by 

xr;(q,c.™) = c/(r/7v) 

Yl Yl *^^j(^ + ^' ^" + Wm)G'j»(-k, -e„) (A7) 

k n 

X^'(q,c.„0 = -C/(T/iV) 

X! Y '-^y (^ + ^' ^" + Wm)Gjj(k, e„). 



(A8) 
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The components of the self-energy matrix are obtained 
from Eq. (0), 
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